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ABSTRACT 

We estimate flux ratios of the extrasolar planet XO-lb to its host star XO-I at 3.6, 4.5, 5.8 and 8.0 
microns with the IRAC on the Spitzer Space Telescope to be 0.00086 ± 0.00007, 0.00122 ± 0.00009, 
0.00261 ± 0.00031 and 0.00210 ± 0.00029, respectively. The fluxes are inconsistent with a canonical 
cloudless model for the thermal emission from a planet and suggest an atmosphere with a thermal 
inversion layer and a possible stratospheric absorber. A newly emerging correlation between the 
presence of a thermal inversion layer in the planetary atmosphere and stellar insolation of the planet 
(jBurrows et al.l[2007b( l is refined. The sub-stellar point flux from the parent star at XO-lb of ~ 0.49 
X 10^ erg cm ~^ s ~^ sets a new lower limit for the occurrence of a thermal inversion in a planetary 
atmosphere. 

Subject headings: stars:individual(XO-l) — binaries: eclipsing — infrared:stars — planetary systems 
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1. INTRODUCTION 

Over 270 extrasolar planets have been reported, 
more than 30 of which transit their primary star^. 
In addition to the mass, radius and inclination of 
the planet evident from transits, atmospheric com- 
position can also be studied through transmission 
spectroscopy, l eading to detections of sodium, water 
and r nethane fCharbonneau et all 120021 : iTinetti et al.l 
120071: ISwain et all |2008[ ). In addition, secondary 
eclipse observations p rovide broadband emission spectra 
(jKnuts on et all 120081 : ICharbonneau et al.' '2008) , plane- 
tary brightness tem p eratures (Charbonncau et al. 200^ 
iDeming et all 120051: [Harrington et al.l I2007D ari d even 
day-night te i npera ture contrast (jKnutson et al.l l2007f ). 
[Torres et al.l ()2008[ ) provide a reanalysis of light curves 
and RV measurements of all then known transiting 
planets. 

For "hot Jupiters", planets with Porb 10 days, a 
favorable planet-star ratio in the IR allows for direct 
detection of the planetary atmosphere by comparing the 
combined flux from the star and the planet during and 
out of secondary eclipse at the superior conjunction. 
The contrast ratio in the mid-lR (1-10 microns) can 
be higher than 10~^ ( c f. ob servations of HD 189733b 
by ICharbo nneau et al] (120081)) a nd theoretical predic - 
tion s bv [Burrows et all (I2006D . iFortnev et all (|2006f ) 
and iBurrows et an ~( 2007al ). "allowing for detection of 
fluxes from planets at secondary eclipse using the 
IRAC, IRS , and MIPS c amera s of the Spitzer Space 
Telescope (jWerner et al.l ()2004[ )). Five planets have 
had their secondary eclipse fluxes measured in one 



or more IRAC bands: TrES-1 



(ICharbonnea u et al.l 
2005[) . HD 209458b (iDeming et al.ll20 05: Knutsonliian 
20081) ■ HD 189733b (IKnutson et al.l " l2007i), HD 149026b 
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(iHarrington et all 12007D and GJ 436b (jPeming et al.l 
|2007| ). In addition low resolution spectra of 2 transit- 
ing planets were obtained with the IR S spectrometer 
betw een ~ 7 and 15 fim: HD 189733b (iGrillmair et all 
I200I and HD 209458b (iRichardson et all 120071 ). 

Recently a detection of an atmospheri c feature at- 
tribut ed to water has be en claimed by iTinetti et al.l 
(|2007| ) and iBarmanI ()2007l ) by studying the transit flux 
ratios of HD 1 89733b and HD 209458b, respectively. 
IBurrows et al.l ()2007a[ ) analyzed the secondary trans- 
mission spectra of HD 209458b at all 4 infrare d IRAC 
Spitzer channels observed by IKnutson et al.l ()2008l ) 
and suggested the observations are consistent with an 
atmospheric thermal inversion layer and yet unknown 
stratospheric absorber. A detailed study of the IR 
secondary eclipse planetary spectra of HD 209458b, HD 
189733b, TrES-1, HD 149 026b and non- e clipsin g HD 
179949b, and v And b by IBurrows et all (|2007bf ) sug- 
gests that the presence of such a stratospheric absorber 
might be dependent on the flux from the star at the 
sub-stellar point on the planet as well as second order 
effe cts like metallic i ty and planetary surface gravity. In 
the IBurrows et al.l ()2007bf ) interpretation planets with 
high sub-stellar point flux (e.g., HD 209458b, OGLE-Tr- 
56b, OGLE-Tr-132b, TrES-2b and X0-3b) would have 
a stratospheric layer and a water feature in emission 
while planets with lower fluxes (XO-lb, TrES-1, X0-2b 
and HD 189733b) wo uld have no such layer and a water 
feature in absorption. iFortnev et al.l ()2007l ) also suggest 
a similar division of planetary spectra based on incident 
stellar flux. Based on th e planetary sub - stellar point 
flux from the star, both IBurrows et al.l ()2007bl ) and 
IFortnev et all (|2007l ) predict that XO-lb should not 
exhibit a thermal inversion in its atmosphere. 

We present observations of the infrared spec- 
tral energy distribution (SE D) of the planet XO- 
lb (|McCullough et all 12001 in all 4 IRAC channels 
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obtained during secondary eclipses with the IRAC of 
Spitzer Space Telescope. By comparing our ^ 4-8 fj,m 
SED with atmospheric models, we test for the presence 
of a thermal inversion layer in XO-lb. 

2. OBSERVATIONS 

The InfraRed Array Camera (IRAC: lFazio et al.ll2004) 
has a field of view of 5.2' x 5.2' in each of its four bands. 
Two adjacent fields are imaged in pairs (3.6 and 5.8 mi- 
crons; 4.5 and 8.0 microns). The detector arrays each 
measure 256 x 256 pixels, with a pixel size of approx- 
imately 1.22" X 1.22". We have observed XO-1 in all 
4 channels in two separate Astronomical Observing Re- 
quests (AORs) in two different sessions: the 4.5 and 8.0 
micron channels for 5.9 hours on UT 2007 April 02 (AOR 
21374464) and the 3.6 and 5.8 micron channels for 5.9 
hours on UT 2007 Aug 10 (AOR 21374208). We used 
the 12 s frame time, obtaining 1620 full-array images in 
each bandpass with a cadence of 13.2 s and an effective 
integration time of 10.4 s . The pointing was not dithered 
and was selected such that for the Apr 2007 observations 
in the 4.5 and 8.0 micron channels, two bright calibra- 
tors for XO-1 (2MASS J16021184-I-2810105 J=9.939 J - 
K =0.412) were used: 2MASS J16021795+2813328: J = 
9.913 J - K = 0.564 and 2MASS J16020133+2809268: 
J = 12.542 J - K = 0.795. The Aug 2007 observations 
in the 3.6 and 5.8 micron channels used 2 bright cali- 
brators: 2MASS J16021311-H2809004: J=11.045 J - K = 
0.652 and 2MASS J16020133-I-2809268: J = 12.542 J - 
K = 0.795. 

We used the standard IRAC Basic Calibrated Data 
products (version 16.1) described in the Spitzer Data 
Handbook^, which includes dark frame subtraction, mul- 
tiplexer bleed correction, detector linearization, and 
flat-fielding of the images. We converted the times 
recorded by the spacecraft in the FITS file header key- 
word DATE-OBS to heliocentric Julian dates using the 
orbital ephemeris of the spacecraft provided by the Hori- 
zons Ephemeris System''. 

Prior to performing aperture photometry, we resam- 
pled the images in all four channels to a 10 times finer 
grid in each spatial direction using flux-conserving b i- 
linear interpolation (similar tol Harrington et"aI1 (|2007f l). 
With the implementation for aperture photometry that 
we used, resampling makes a marginal improvement in 
the photometry, i.e. a slightly lower r.m.s. of out-of- 
eclipse points, presumably related to how the routine 
handles fractional pixels at the edge of the aperture. 

The zodiacal background was subtracted in each 
channel by constructing a histogram of all pixels in each 
image and fitting a Gaussian to the distribution of the 
zodiacal background brightness. A constant mean value 
of the Gaussian was then subtracted from each pixel in 
the image to construct a background-subtracted image. 
The centroids of XO-1 and the 2 calibrators were evalu- 
ated by fitting a Gaussian to the stellar flux distribution. 
The pointing varied by 0.3 pixel, and the shifting of the 
stellar centroid within a pixel, which have sub-pixel sen- 
sitivity variations, resulted in a modulation of the stellar 
flux in the 3.6 and 4.5 micron channels (described below). 

® http: //ssc ■ spitzer ■ caltech. edu/lrac/dh/ I 
^ http : //ssd ■ jpl ■ nasa . gov/ | 



Aperture photometry was then performed on the im- 
ages with an aperture radius of 4 pixels, which was found 
to be the optimum value for all 4 channels. The size 
of the aperture was determined by minimizing the rms 
scatter in the light curve for observations outside of the 
eclipse. Apertures smaller than 4 pixels contained insuf- 
ficient stellar flux and larger apertures were more con- 
taminated by the sky (especially in the high background 
signal of the 5.8 and 8.0 micron channels). An appro- 
priate aperture correction for each channel was applied 
to the stellar flux value according to the Spitzer Data 
Handbook of [1.112, 1.113, 1.125, 1.218] for the [3.6; 4.5; 
5.8 and 8.0] micron channels, respectively. The a of the 
out-of-eclipse points was calculated iteratively using 3-tT 
outlier rejection at each step until no more points were 
rejected. To remove cosmic rays the resultant robust a 
was used to reject entire images which contain the 3-ct 
outliers above and below the mean of the light curve. 
1.4%; 1.9%; 3.0 %, 2.4 % of images from the 3.6, 4.5, 
5.8 and 8.0 micron channels, respectively, were removed 
in this fashion. The higher rejection rate in the two red- 
der channels is consistent with a high er number of cosmi c 
ray-affected pixels in these channels (jPatten et al.ll20d^ . 
Throughout the analysis we have preserved fiux units. 

2.1. 3.6 and 4-5 micron time series 



The 3.6-/im time series exhibited a sharp increase dur- 
ing the first ^--^30 minutes of exposure for XO-1 and the 
2 calibrators, presumably as a result of the instrument 
reaching a new equilibrium after previous observations. 
Such relaxation effects can reach several percent and usu- 
ally stabilize within the first hour of observations of a new 
target. We have ignored the first 125 points (~30 min) 
in the 3.6-/im time series in addition to the high-sigma 
outfier rejection as described in the previous section. 

A strong correlation between the sub-pixel centroid 
and stellar brightness was observed in both the 3.6- 
fim and 4.5-/im channels, with flux magnitudes of ^ 
0.6% and 0.8%, respectively. This well studied effect 
(|Charbonneau et al.l[2005l : lMorales-Calder6n et al.ll2006l ) 
is due to the InSb detector intrapixel sensitivity vari- 
ations as the spacecraft jitters ^--^0.3 arcsec in orienta- 
tion over a period of ~3000 seconds ^. The uncorrected 
sub-pixel intensity variations are clearly visible in the 
time series of XO-1 in the 3.6-/im and 4.5-/im channels 
in Fig. [TJ We h ave corrected for thi s sub-p ixel intensity 
variations after ICharbonneau et"aI1 (|2008f l by fitting a 
quadratic function to the photometric flux points of XO- 
lb observed out-of-eclipse as a function of the x and y 
sub-pixel centroids: 

Isubpixei = 6i+ 62 X 2; + 63X biX y + b5 X y^, (1) 

where x and y are the subpixel centroids of center of 
light of the star and hn are the fit parameters. The rms 
residual of the XO-1 time series for points outside of the 
eclipse after correction for the sub-pixel intensity vari- 
ation was 0.0020, which is 18% higher than a theoret- 
ical estimate of XO-1 Poisson noise based on detector 
read noise and background noise. The time series was 



: //ssc . spitzer . caltech . edu/document s/exoplanetmemo . txt 
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TABLE 1 

XO-1 ABSOLUTE FLUXES 



-0.10 -0.05 -0.00 0.05 0.10 
Time from predicted eclipse center (days) 



0.15 



Fig. 1. — Secondary eclipse observations of XO-lb with 
IRAC on Spitzer Space Telescope obtained on UT 2007 April 
02 and UT 2007 Aug 10 in 3.6, 4.5, 5.8 and 8.0 micron chan- 
nels (from top to bottom) binned in 6-minute interval and 
normalized to 1 and offset by 0.02 for clarity. The overplot- 
ted solid lines do not represent a fit to the data, but rather 
show the correction for the detector effects. The 3.6 and 4.6 
micron time series are decorelated using XO-1 out-of-eclipse 
points and the 5.8 and 8.0 micron time series is detrended 
using a fit to a calibrator star in the field (see i|2.1l fc i|2.2l for 
details.) 

then normalized by taking the robust average of out-of- 
eclipse points and binned in 6-minute intervals contain- 
ing approximately 30 individual measurements each (see 
Fig. El). 

The 4.5-/im time series also exhibited an initial 
relaxation-induced brightness increase and consequently 
139 points corresponding to the first ^--^30 minutes of ob- 
servations were rejected, which is more points than in 
the 3.6-^m time series. The analysis of the time series 
was identical to that of the 3.6-/im time series. The rms 
of out-of-eclipse points was 0.0024, which is 19% higher 
than the theoretical estimate and is similar to that of 
TrES-1 (rms=0.0027 lCharbonneau et all (l2005l )). 



IRAC channel 
effective A 
(microns) 



XO-1 ffux 



(mjy) 



XO-1 instrumental 
magnitude 
(mag) 



3.6 
4.5 
5.8 
8.0 



45.1 ± 2.3 
28.5 ± 1.4 
18.1 ± 0.9 
11.7 ± 0.6 



9.49 
9.50 
9.51 
9.35 



We have tested tor tiie Imearity of the detector re- 
sponse in the 3.6 and 4.5 micron channels in which 
XO-1 is close to the onset of detector non-linear re- 
sponse. Using a sub set of data from the SAGE survey 
(|Meixner et al.|[2"006D obtained in the high dynamic range 
(HDR) mode of IRAC camera with both 0.6s and 12.0s 
integration times we are able to determine that both the 
3.6 micron and 4.5 micron XO-1 fluxes are unsaturated 
and in the detector linear regime response. Table[l]shows 
the absolute XO-1 fluxes and the instrumental magni- 
tudes in the four IRAC channels. 

2.2. 5.8 and 8.0 micron time series 

The 5.8 and 8.0 micron time series were recorded with 
Si:As detectors and do not thus exhibit the prominent 
sub-pixel intensity variations evident in the 3.6 micron 
and 4.5 micron channels. The first ~30 minutes of ob- 
servations (139 data points) were rejected as the instru- 
ment settled into a new equilibrium state. Fig. [1] shows 
intensity variation with time, which is caused by changes 
in the effective gain of individual pixel s over time. This 
effect has also been been observed by iHarrington et al ' 
(2007') at 8 micron with IRAC and bv iDeming et al 



2006) at 16 micron with IRS. The intensity variations 
are depen dent on the illumina t ion le vel of the individ- 
ual pixel (jKnutson et al.l l2007l l2008f ). pixels with high 
illumination will reach their equilibrium within ~1 hour, 
but lower illumination pixels increase in intensity over 
time, approximately proportional to the inverse of the 
logarithm of illumination. We have decided not t o cor- 
rect each pixel in the image as lKnutson et al.l ()2007f ) have 
done in their 33-hour observation of HD 189733 in the 
8 micron channel of IRAC and instead fit a combined 
linear and quadratic logarithm function of time from the 
beginning of observations to the time series of the bright 
calibrator 2MASS 16020133+2809268. The fit to the 
time series of the calibrator was then used to remove the 
detector ramp from the time series of XO-1 as depicted 
in Fig. [TJ 

Imodei = Cii+a2'x Ai -|- 03 X InAt + a^x (InAt)'^, (2) 

where Imodei is the model fiux. At is the time since the 
beginning of observations and are the free parameters. 
The detector ramp intensity decreased in fiux during the 
5 hours of ob servation by ~ 0.2 % , following a similar 
trend seen bv iKnutson et all (|2008l ) in their 5.8 micron 
time series of brighter HD 209458b. After removal of the 
detector ramp and normalization, the rms of unbinned 
5.8 micron out-of-eclipse points for XO-1 was 0.0081, 
which is 44% higher than theoretical Poisson noise, based 
on detector read noise and background noise. ^ The nor- 

^ Ihttp : //ssc. spitzer . caltech . edu/document s/irac_memo . txt I 
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malized and binned 5.8 micron time series is depicted in 
Fig.H 

Fitting the detector ramp using a calibrator star, which 
is non-variable at the 0.012 level after detector ramp re- 
moval, allows us to bypass using XO-1 itself to remove 
the detector ramp by making a fit to its 5.8 micron out- 
of-eclipse points. The choice whether to correct for the 
detector ramp of XO-1 using the photometry of a calibra- 
tor star or the target star itself could be a limiting factor 
in our analysis. Deming (2008, p.c.) has indicated that 
in a long (> 10 hours) series of full-frame 12s-exposure 
IRAC photometry, fitting the detector ramp to the tar- 
get star while masking out the points in the eclipse could 
be more appropriate. We implemented each of the two 
alternate methods of correction, and for these observa- 
tions of XO-lb, the measured depths at 5.8 micron are 
w ithin 1-a of each oth er. 

ICharbonneau et all (I2005D [iTnutson et all (I2007D . 



iHarrington et all (|2007f) and iDeming et al.l (|2007h have 
reported a nonlinear flux increase over time in the 8.0 mi- 
cron IRAC channel. We also detect a nonlinear increase 
in the brightness of XO-1 and in the 2 calibrators, in 
addition to the sharp increase during the first ~30 min- 
utes of observations (137 points). The initial ramp-up 
data were discarded and the ^ 1.6% non- linear increase 
of the XO-1 time series over 5 hours (see Fig. [T]) was 
removed by fitting a combined linear and quadratic log- 
arithm function (eq. [5]) to the time series of the first 
calibrator 2MASS J16021795-f2813328 and dividing the 
XO-1 time series by the fit. The time series was normal- 
ized using out-of-eclipse points and binned into 6-minute 
bins (Fig. [2]) for viewing clarity. The resulting rms of 
unbinned out-of-eclipse points for the XO-1 time series 
was 0.0075, 74% above Poi sson noise, but similar to that 
of TrES-1 (rms = 0.0085, (|Charbonneau et al.l[2005h '). 

3. ANALYSIS 

We fit th e secondary eclipse ligh t curves using the for- 
malism of iMandel fc Ago! (|2002[ ) with no stellar limb 
darkening and adop t stellar and orbital parameters 



5+0.018 



R* — 0.928_o;oi3 R-O: R-p 

I +0.46 



dHohnan e Fall [2001 ^°: 

1.184;g°28 R^^pii^ i ^ 89.3i;^;5;^ degrees, and a 
0.0488 ± 0.0005 AU with ephemeris ()McCullough et all 
[2006.1 



TciE) = 2, 453, 808.9170(ff JD) + i;(3.941534 days) . 

(3) 

We fit the depth of the eclipse AF and the timing of 
the centroid AT independently in all 4 channels in the 
unbinned li ght series using Levenberg-Marquardt mini- 
mization (Press et al.lll9 92) with an equal error assigned 
to all points, which is equal to the rms of out-of-eclipse 
points in each time series. Best-fit eclipse curves are 
plotted in Fig. [2] and the eclipse parameters are listed in 
Tabled] They are the channel wavelength, eclipse depth 
AF, eclipse mid-center time in HJD, and the timing off- 
set At in minutes from the expected secondary eclipse 
mid-center time for an assumed eccentricity of zero, and 

We have also reduced the dat a using the stellar and orbital pa- 
rameters from [Torres et all l|2008l ) as a test, but the eclipse depths 
changed negligibly and eclipse mid-center timings were all within 
1-a. 

" 1 Rjup = 71,492 km. 
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Fig. 2. — Secondary eclipse of XO-lb observed with IRAC on 
Spitzer Space Telescope in 3.6, 4.5, 5.8, and 8.0 micron chan- 
nels (top to bottom) corrected for detector effects, normalized and 
binned in 6-minute intervals and offset for clarity. The best-fit 
eclipse curves are overplotted. 

the reduced ■ The reduced is close to 1.0 in all 4 
channels, indicating a good fit to the data. 

To estimate the errors on the depth and mid-eclipse 
timing we performe d the error a nalys is using the boot- 
strap method from iPress et al.l ()1992f ). The bootstrap 
method makes no prior assumptions about the distribu- 
tion of the noise in the data and the data points are not 
altered as in the Monte-Carlo analysis. For 10,000 trial 
runs we have randomly drawn with replacement points 
from the normalized- and detector-effect-decorrelated, 
but otherwise unaltered, light curve, until we had the 
same number of data points in the light curve that 
we started with. During each iteration we performed 
the full eclipse fitting for eclipse depth AF and eclipse 
mid-center AT. The 1—a errors for AF and AT were 
computed by fitting a Gaussian to the respective 1-D dis- 
tribution of bootstrap points and are reported in Table[21 

The eclipse depth errors AF = [0.00007, 0.00009, 
0.00031, 0.00029] for the J = 9.939 XO-1 compare fa- 
vorably with the eclipse depth errors [0.00009, 0.00015, 
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TABLE 2 

Secondary eclipse best fit parameters 



A 


Eclipse Dcptli AF 


Eclipse Center Time 


Time offset AT 


Reduced 


(microns) 




(HJD) 


(min) 




3.6 


0.00086 ± 0.00007 


2454323.28342 ± 0.00242 


-5.42 ± 3.49 


0.76 


4.5 


0.00122 ± 0.00009 


2454193.21175 ± 0.00259 


-6.94 ± 3.72 


0.84 


5.8 


0.00261 ± 0.00031 


2454323.29066 ± 0.00488 


5.00 ± 7.03 


0.89 


8.0 


0.00210 ± 0.00029 


2454193.22043 ± 0.00486 


5.56 ± 7.00 


0.99 



0.00043, 0.00026] of iKmitson eTall ([2008h for the J = 
6.591 HD 209458, despite the fact that XO-1 is dimmer. 
XO-1 observations were made in the fuU-array mode 
with 10.4-s integration time and readout time 2.8 s for a 
total of 1,620 images, while the HD 209458 observations 
were made in sub-array mode with exposure time of 0.1 
s in sets of 64 in each channel totaling 26,880 usable 
images. The S/N scales as cx y/n^ * fe * At / atotai, 
where Uexp is the number of exposures during the 
duration of the eclipse, /e is the stellar flux in signal 
electrons, At is the integration time and atotai is the 
combined Poisson, readout and background noise. The 
predicted eclipse depth errors for XO-1 in the four IRAC 
channels are then [0.4; 0.5; 0.9; 2.5] times the respective 
eclipse depth errors for the 21 times brighter HD 209458. 

To test the robustness of our data reduction and 
analysis technique and consistency with other observa- 
tions in the IRAC full-array mode we have re-reduced 
the 4.5 and 8.0 micron IRAC seco n dary eclipse data 
of TrES-1 by ICharbonneau et all (|2005D with our 
pipeline, rejecting the first 30 minutes in both channels. 
The procedure described in ^2.11 and ^2.21 was used 
to gether with up dated stellar and planetary parameters 
of [Torres et al.l (|2008) to derive the eclipse depths 
relative to the star TrES-1 of AF4.5;.m = 0.00043 and 
AFs.Ofj,m = 0.00194, which a r e -l.l g and -0.9cr away 
from the ICharbonneau et"aI1 (|2005f ) secondary eclipse 
depths of AF4.5;.m = 0.00057 ± 0.00013 and AFg.o^m 
= 0.00225 ± 0.00036, respectively. Our secondary 
eclipse [4.5; 8.0] micron channel mid-center timing offsets 
of [ +12.5; -2.3] minute s are [-l.lg; -2. Oct] away from 
the ICharbonneau et al.l (|200l [4.5; 8.0] micron channel 
values -1-19.6 ± 6.6 min and -1-8.3 ± 5.2 minutes, 
respectively. While the eclipse depths are consistent at 
the ~ 1-CT level, the eclipse mid-center timing offset in 
the 8.0 micron channel is on ly mildly consistent at the 
2. Oct level, probably because ICharbonneau et al] ()2005f ) 
observations were made using multiple Astronomical 
Observation Requests (AORs) whi ch resulted in arbi- 
trary flux shifts in the time series. ICharbonneau et al.l 
(|2005f ) do not mention how they have corrected for these 
flux shifts. Despite this fact our pipeline is capable 
of reproducing their results to within ^ 1-ct in eclipse 
depth and ~ 2-ct in mid-eclipse timing. We have thus 
demonstrated that our reduction pipeline is robust and 
the secondary eclipse depth estimates are consistent 
with other major full-array pipelines. 



4. DISCUSSION 

The eclipse mid-center timings for XO-lb in Table 
[2] are individually consistent with zero timing resid- 



uals for a circular orb it based on the ephemeris by 
iMcCullough et al.l (|2006[ ) . but the April and August 2007 
combined timings show a time shift. The UT Apr 2 2007 
observations of 4.5 micron and 8.0 micron channels have 
a combined eclipse mid-center timing of 2454193.21366 ± 
0.00228 HJD, with a delay of -4.2 ± 3.3 min, while 
the UT 2007 Aug 10 observations of 3.6 micron and 
5.8 micron channels have a combined center of eclipse 
time of 2454323.28485 ± 0.00217 HJD, with a delay of 
-3.4 ± 3.1 min. The mid-eclipse timing offsets can be 
interpreted as evidence for further planets in the sys- 
tem, as can a non-uniform brightness distribution of the 
planet (which can change the shape of eclipse ingress 
and egress and, thus , shift the eclipse mid-center time 
(|Williams et al.ll2"006D ). or as evidence for a non-zero ec- 
centricity. Using the larger mid -eclipse April time offset, 
the approximate equation (e.g. iKopall (|1959f) Eq. 9-23) 



e X cos{uj) 



-KAt 



(4) 



where e is the eccentricity, uj is the longitude of perias- 
tron, P is the orbital period, and At is the centroid time 
shift from expected time of secondary eclipse, allows us 
to set a 2 CT upper limit on excos{Ljj) < 0.0024. 

The XO-lb eclipse depths in Fig. [3] show several 
trends. The planet-to-star contrast of XO-lb peaks in 
the 5.8 fim channel, with a decrease in the 3.6 micron 
and 4.5 micron channel and a slight decrease towards 
the 8.0 micron channel. Furthermore the flux in the 4.5 
fim channel is higher than in the 3.6 micron channel, 
which does not match the gene r al cha racter of the cloud- 
less models of iBurrows et al.l (|2006t ) for redistribution 
parameter P„ =0.3 (dot-dashed hne and open circles 
for band-averaged ratios in Fig. |3|), which predict a 
lower flux at 4.5 microns than at 3.6 microns; the model 
with thermal inversion predicts the opposite. Since the 
observations manifest a higher flux in the 4.5 micron 
channel than at 3.5 micron, a thermal inversion in the 
atmosphere might be indicated. Pn=0 corresponds to 
no heat redistribution from the planetary day-side to 
the night-side and Pp = 0.5 stands for full redistribution 
(see IBurrows et al"] ()2007b[ ) for details). The possibihty 
of thermal inversi on in a plan e tary a tmosphere has 
been suggested by 'Hubenv et al.l ( 2003fl. IBurrows et al.l 
2006, 2007a) and Fortne v et al l (1200711 Recently, 



Burrows et al.l (12007a. b) suggested that model spectra 
could match the observations of HD 209458b and 
HD 149026b if a stratospheric absorber of unknown 
composition (possibly tholins, polyacetylenes, TiO or 
VO) were present in the atmosphere of the planet. 
The presence of a stratospheric absorber would yield a 
thermal inversion in the planetary atmosphere and the 
presence of the water features in emission for a variety 
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of heat redistribution parameters P„ 




Wavelength tim 



Fig. 3. — Spitzer Space Telescope IRAC secondary eclipse 
depths for XO-lb with bootstrap error bars (filled squares). The 
predicted emission spectrum of the planet with an upper atmo- 
spheric absorber of Ke = 0.1 cm^/g and a redistribution parameter 
of P„=0.3 is plotted as a solid line. A model with no atmospheric 
absorber and a redistribution parameter of P„=0.3 is over plotted 
with dot-dashed line (see [[3] for details). The band-averaged flux 
ratios are plotted as open squares and open circles for the models 
with and without stratospheric absorber, respectively. The normal- 
ized Spitzer Space Telescope IRAC response curves for the 3.6-, 4.5- 
, 5.8-, and 8.0 micron channels are plotted at the bottom of the fig- 
ure (doted lines). The theoretical flux ratios obtained from a XO-1 
stellar spectrum (from http : //kurucz .harvard. edu/stars/XOl/ I 
and an assumed black-body spectrum tor the planet at [12U0, IfciOU] 
K are plotted as dashed lines. 



Our observations suggest the presence of a thermal 
inversion layer and a possible stratospheric absorber 
in the atmosphere of the XO-lb planet. The solid line 
and open squares in Fig. [3] depict an atm ospheric model 
of XO - lb, following the methodology of iBurrows et al.l 
()2006l 12007a ) with a thermal inversion and a strato- 
spheric absorber of opacity of — 0.1 cm^/g and 
redistribution parameter of P„ — 0.3. The latter model 
fits the data better than the canonical cloudless model 
with P„ = 0.3 (dot-dashed curve and open circles for 
averaged band ratios). The band-averaged flux ratios 
for the model with a stratospheric absorber (open 
squares) are within the error bars for the 3.6, 4.5, and 
8.0 micron channels, but are inconsistent by 2.7a with 
the band-averaged flux ratios for the 5.8 micron channel. 
This is similar to the situation f or the IRAC fit to th e 
observations by HD 209458b bv iKnutson et all (|2008D . 
The absorber-free canonical model (dot-dashed line) 
is clearly inconsistent with our observations (Fig. [3]) 
of XO-lb in all 4 channels by [3.4cr, 7. la, 6.3a, 3.7a], 
respectively. 



IBurrows et all (|2007b( l and iFortnev et all (|2007| ) sug- 
gested that the presence of the stratospheric absorber 
might be correlated with the incident flux from the 
star at the sub-stellar point on the planet, the precise 
level of which is yet to be refined. The presence of 
an irradiation - induced strato s phe ric absorber has been 
suggested by IBurrows et all JlOOTa) for HD 209458b 
(see our Fig. |4|) with a sub-stellar flux of ~1.07 x 10^ 
erg cm ~^ s at a distance a — 0.045 AU. Interestingly, 
XO-lb has a lower sub-stellar flux of ~ 0.49 x 10^ erg 



cm s and a semi-major axis of a = 0.0488 AU, 
but still manifests evidence for a thermal inversion. A 
recent study of the broadband in frared spectrum of 
HP 1 89733b (see our Fig. \^ by ICharbonneau et all 
(2008) finds no evidence for an atmospheric thermal 
inversion, despite a simila r sub-stellar point flu x of ~ 
0.47 X 10^ erg cm s ([Burrows et al.ll2007bD with a 
smaller semi-major axis a = 0.0313 AU. Further study 
of planetary atmospheres should refine the concept of 
this sub-stellar flux boundary with respect to the pres- 
ence/absence of a stratospheric absorber and thermal 
inversion. 



: .XO-lb 

0040 - • HD 209458b 

- 4 HD 189733b 

- ♦GJ436b 
: T TrES-1 

- ,HD 149026b 



.A 




Wavelength 

Fig. 4. — Comparison of Spitzer Space Telescope IRAC secondary 
eclipse depths: XO-lb (filled square) from this paper; HD 209458b 
(filled circle ) from Knutson ct al. (2008); HD 189733b (filled tri- 
angle) from JCharbonn cau ct al. (2008); GJ 436b (filled star) from 
Deming et al.l II2007I) ; TrES-1 (filled upside down triangle) from 
Charbonneau et al.l ||200 5); and HD 149026b (open circle) from 
Hsgrmgton'er^l.l 1120071 ). The central wavelengths have been offset 
by [-|-0.1;-0.1] microns for clarity. The normalized Spitzer Space 
Telescope IRAC response curves for the 3.6-, 4.5-, 5.8-, and 8.0 mi- 
cron channels are plotted at the bottom of the figure (doted lines). 



Atmospheric water detection has been claimed in 
the transit broadb and spectra of HD 189733b by 
iTinetti et al.l (I2007D and in its secondary eclip s e spec- 
tra bv IFortnev fc MarlevI (|2007D . IBurrows et all (|2007af) 
also found evidence for water vapor emission in the at- 
mosphere of HD 209458b. Our data can be interpreted as 
evidence for rovibrational band of water emission long- 
ward of ~4.0 microns, which manifests itself as a flux 
enhancement in Fig. |3| compared to the cloudless model. 
The depth of the flux ratio trough near the 3.6 micron 
channel in the atmospheric models with stratospheric ab- 
sor ber strongly depends on the redistribution parameter 
P„ (jBurrows et al.l ()2007bl ). especially their Fig. 4). Fur- 
ther modeling would allow tighter constraints on the P„ , 
not just for XO-lb, but for a variety of planets. The 
fortuitous importance of the 3.6 micron IRAC channel 
to the study of planetary atmospheres is likely to be en- 
hanced as the Spitzer Space Telescope runs out of cryo- 
coolant in 2009, when only the 3.6 micron and 4.5 micron 
channels will be available. 

5. CONCLUSION 

We report the estimated flux ratios of the planet XO- 
lb in the Spitzer Space Telescope IRAC 3.6, 4.5, 5.8 and 
8.0-fim channels. We find that the estimated fiuxes are 



XO-lb Thermal emission 



7 



not consistent with a canonical cloudless model for ther- 
mal emission from the planet and instead may indicate 
an atmosphere with as yet unknown stratospheric ab- 
sorber and a likely thermal inversion, which would cause 
the water band longward of 4.0 microns to switch from 
absorption to emission. The atmospheric model with a 
thermal inversion produces a tight match to the data at 
3.6, 4.5, and 8.0 microns, but is inconsistent by 2.7a with 
observations at 5.8 m icrons. This is sim ilar to observa- 
tions of HD 209458b (jKnutson et al.1 '20081. 

The presence or absence of the stratospheric absorber 
and thermal inversion layer has been linked to the flux 
from the parent star at the sub-stellar point on the 
planet. The XO-lb sub-stellar point flux of ~0.49 x 
10^ erg cm ~^ s ~^ is the lowest so far reported for a 
planetary atmosphere with a thermal inversion. Obser- 
vations of atmospheres of other planets may permit a 
better understanding of the thermal inversion layer and 
parametrization of the characteristics that create such a 
thermal inversion. 
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